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The small oscillations of solitons in 2D Bose-Einstein condensates are inves- 
tigated by solving the Kadomtsev-Petviashvili equation which is valid when 
the velocity of the soliton approaches the speed of sound. We show that the 
soliton is stable and that the lowest excited states obey the same dispersion 
law as the one of the stable branch of excitations of a ID gray soliton in a 
2D condensate. The role of these states in thermodynamics is discussed. 



Solitons in Bose-Einstein condensates (BECs) have recently attracted 
much attention and have been studied both theoretically and experimentally. 
The Gross-Pitaevskii (GP) equation admits solitonic solutions corresponding 
to a local density depletion (i.e., gray and dark solitons). They have been 
created in 3D BECs 1 and their decay into vortex rings has been observecP. 
In 2D, an interesting type of soliton corresponds to a self-propelled vortex- 
antivortex pair that, when it moves at a velocity close to the Bogoliubov 
sound speed, c, takes the form of a rarefaction pulsed The energy E and 
momentum P were calculated in RefP, and it was found that the curve E(P) 
is below the Bogoliubov sound line, approaching it in the low momentum 
limit. We investigate the excitation spectrum. Here we present our first 
results for solitons moving with velocity close to c, for which one can reduce 
the GP equation to the simpler Kadomtsev-Petviashvili (KP) equatiori^. 

Let us consider a 2D condensate with a soliton moving at velocity V in 
the x-direction. If the density at large distances is n^, one can define the 
healing length £ = h/(2mgn 00 ) 1 ^ 2 , where g is the mean-field coupling con- 
stant and m is the mass of the bosons. One can introduce the dimensionless 
variables x — > £x, y — > £y and t — > mt^ 2 /K, the normalized order parameter 
* — ► \J n oo^ and the velocity U = m^V/fi = V/(cy/2) where c = h/(y/2m£) 
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is the sound speed. In the frame moving with the soliton, the GP equation 
becomes^ 

2i— * = _V 2 ^ + 2i[/— tf + (|^| 2 - l) . (1) 

The order parameter can be written as ^ = n 1 / 2 e* 5 . In the limit of V — > c, 
one can introduce the small parameter e = y/l — 2U 2 , so that U — l/y^2 — 
l/(2\^2)e 2 . The density depletion associated with the soliton becomes shal- 
low and one can expand the density and the phase as 

n = l-e 2 f + ... (2) 
S = es + ... . (3) 

To the lowest order in e, the GP equation gives d x s = — ef/y/2 and the 
function / obeys the KP equation 

d fdf | df | d 3 f\ _ d 2 f 
dx I dt dx dx 3 J dy 2 ' 

where we have introduced the stretched variables x = —ex + £ 3 t/(2\/2), 
y = e 2 y/V2 and t = e 3 t/(4 v / 2). 

i) ID gray soliton in 2D condensate. 

By assuming the stationary solution of (J3J to be independent of y, 
Eq. (@J) reduces to the Korteweg-de Vries (KdV) equation and the solu- 
tion is fo(x — 2t) = sech 2 [(x — 2t)/y/2], corresponding to a gray soliton. 
The linear stability of this solution can be studied by taking f(x, y, t) = 
fo(x — 2t) + ip(x — 2t)e l ( k y~ ult ' > and linearizing Eq. (@J) in ip{x). One gets 

d 4 ip d 2 , s d 2 4> ~ 9 dip 

This equation was already solved by Zakharo^, by applying the inverse scat- 
tering method, and also by Alexander et a/P. The excitations which are lo- 
calized along x and periodic along y have dispersion law Co 2 = (16Ar/3\/3)(A;— 
v / 3/2). For a given e, the eigenfrequencies u in the original units of the GP 
equation Q are given by u = £ 3 cD/(4\/2), and the wavevector transforms as 
k = e 2 k/V2, so that 

;„2 / rs" \ 

(6) 

For k > e 2 ^/3/8 the frequency is real, while for k < e 2 ^/3/8 it is imaginary, 
causing the instability of the soliton against bendingPln the limit e — > 0, the 
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Fig. 1. Profile of the 2D soliton, Eq. (JJJ), (upper left) and three examples of 
excited states with Q = 5.0 (upper right), Cj = 11.8 (lower left), and Cj = 24.4 
(lower right). For the excitations, the real part of the odd eigenfunctions is 
shown. 

region of instability becomes vanishingly small and, for k e 2 , the stable 
branch of excitations becomes u = (2/27) 1 / 4 A; 3 / 2 . It is worth noticing that 
this dispersion exhibits the same power law as for capillary waves on the 
surface of a liquid. 



ii) 2D soliton. 

The solution of Eq. Q corresponding to a 2D soliton ia^ 



f (x-2t,y)=4 



+ 2y 2 



(x 



2tf 



[3 + 2 y2 + ( £ _ 2t)2]2 



(7) 



This function is plotted in the upper left frame of Fig. 1. In the units of the 
GP equation l|T|). the density profile of the soliton is found by inserting ( |7| ) 
into PJ). Its size along x and y is of the order of 1/e and 1/e 2 , respectively. 
To the first order in e, the momentum P and the energy E per unit length 



ar 



e 3 



8/T 



8tt 



P = —V^hrioc^e , E 1 = — h 2 nooe = cP 



3m 



(8) 



As it must be, the soliton has a sound-like dispersion in the limit where KP 
equation is valid, that is, when P — > 0. 
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Now we consider fluctuations of the form f(x,y,t) = fo(x — 2t,y) + 
tp(x — 2t,y)e~ tujt and linearize Eq.© by ip. Thus, ip(x,y) satisfies 

— + 6— ff -0)-2— - — = iw— (9) 
dx A dx 2 dx 2 dy 2 dx 

Looking for bound excited states (i.e., ip — > at large distances), it is con- 
venient to integrate Eq. @ twice along x and introduce the function <p(x, y) 
such that ip = The equation for <fi is 

<*± + 6/ ^ - 2^ - ^ = iu^ (10) 
dx 4 dx 2 dx 2 dy 2 dx 

We numerically solve this equation in a square box of size L, by expanding (p 
as 4>{x,y) = J2^4>u^Xi,u(x)x2Ay) with XvO 2 ) = L~ 1 / 2 e j27rra / L for |z/| < I 
Concerning X2,^(y) we note that Eq. ® is invariant for y — > — y, so that the 
function (j) is either even or odd function of y. If it is even, then one can 
take X 2,M = (VL) 1/2 cos(2vr M y/L) for 1 < fjt < I and X 2,o(v) = (V^) 1/2 - 
If it is odd, one can take xzAv) = ( 2 / L ) 1/2 sin(27r//y/L) for 1 < p < I. One 
thus obtain the following matrix equation 

(~Qx ~ 2 Qx ~ Qy/Qx) 4>utx + 6 ^x M vti,v>ti'^v'n' = &4> vtx , (11) 

v',n> 

with q x = 2irv/L, q y = 2-Kfi/L and 

[L/2 

M^y^' = / dxdy Xi,uX2,nfoXiyX2y ■ (12) 

J — Z//2 

The size of the matrix is N x N, where N = 21(1 + 1) for 4> even, and N = 2l 2 
for 4> odd. Typical values in our calculations are I = 70 and L = 60. 

It is found that all the eigenvalues are real. This is consistent with the 
results of RefsP^ and shows the linear stability of the 2D soliton. One can 
expect this result also from pure kinematic considerations. An excited state 
is unstable if it can decay into several phonons. Energy and momentum 
must be conserved in this process. It was proved by Iordanskii and one of 
the authorPI that for a Bogoliubov-like spectrum the conservation laws can 
be satisfied only if the dispersion law E(P) is above the sound line. Since 
the dispersion law of the 2D soliton is below the sound line, the soliton is 
stable. 

We find stable localized modes with positive eigenvalues. Examples are 
shown in Fig. 1 for lD = 5.0, 11.8 and 24.4. The function ip oscillates along 
x and y, with pronounced maxima in the soliton region. Since the density 
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of a ID gray soliton, Eq. ©. 



profile of the soliton decays algebraically, ip(x,y) decays rather slowly. It is 
found that the real and imaginary part of ift are either even or odd func- 
tions of x and each eigenvalue is degenerate, corresponding to even and odd 
functions of y. 

By looking at the oscillations of the eigenvectors ip in the soliton region, 
one can estimate the transverse and axial wavevectors, k x and k y respectively. 
In Fig. 2 we plot the calculated eigenfrequencies Co as a function of k y for 
the lowest bound states (points with error bars, where the error bars are of 
the order of the inverse of the axial size of the soliton). In the same figure 
we plot the dispersion law © of the stable branch of excitations of a ID 
gray soliton. The spectrum is very similar. The 2D soliton have a discrete 
spectrum of bound states, due to its finite size, and all wavevectors are above 
the threshold for the instability of the ID soliton. 

Notice that the soliton oscillations can be subject to damping by emis- 
sion of sound waves. This means that the energy levels of Fig. 2 have to be 
considered as resonance states.^ However, there are reasons to believe that 
this damping is small and can be neglected both in the calculation of the 
energy levels and in the thermodynamic considerations below. We plan to 
investigate this damping in the future. 

As it was noticed in RefP, this 2D soliton can contribute to the specific 
heat C of the 2D Bose gas. Indeed, if one applied the usual Bose statistics 
to the soliton branch, due to its sound- like dispersion (|5|>. its contribution 
would be C oc T 2 as for phonons. The presence of the excited states of the 
soliton can change the situation. The inclusion of these excitations is an 
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open problem; here we just want to provide an argument which illustrates 
its nontriviality. 

One can take into account the excitations of solitons by assuming that 
the solitonic branch of the spectrum depends on two quantum numbers: P 
and k. For example, within the limits of applicability of the dispersion law 
©, one can write 

E(P,k) = c\P\ + (3k 3/2 . (13) 

One has to accept the hypothesis that the number of states in the interval 
d 2 P is, as usual, Sd 2 P/ (2irfi) 2 , where S is the sample area. On the contrary, 
the quantization rule for k is fixed by the length of the soliton in the y 
direction, R y , which behaves as R y oc £/e 2 . Correspondingly, the number 
of states in the interval dk is oc Rydk/(2ir) 3 oc ((,/e 2 )dk. The energy of 
the gas is thus £ oc / d 2 Pdk e- 2 E(P,k)[exp{E(P,k)/T) - l]" 1 . A simple 
calculation shows that the relevant values are E oc e oc T, Pj oc T with 
i = x,y, and k oc T 2 / 3 . This gives the specific heat C oc T 2 / 3 , which 
exceeds the phonon contribution at low T. Of course, this result cannot be 
considered as a rigorous one, in particular because solitons occupy an area 
oc 1/e 3 and they can overlap at low temperature. Nevertheless it shows that 
the problem of the low-temperature specific heat of the 2D Bose gas deserves 
both theoretical and experimental investigation. A direct diagonalization of 
the Hamiltonian or Monte Carlo simulations at finite temperature could be 
suitable tools to test our semi-quantitative prediction. 

We thank C. Lobo, P. Pedri, N. Prokof'ev and D. Thouless for useful 
discussions. 
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